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A theoretical framework is presented, which allows to explain many experi- 
mental facts related to pinning and cross-flow effects between flux tubes in 
type-II superconductors. It is shown that critical state principles, in the 
manner introduced by C. P. Bean for parallel vortex lattices, may be used to 
describe the observed behavior. We formulate a least action principle, giving 
place to a variational interpretation of the critical state. The coarse-grained 
electrodynamic response of the superconductor is solved by minimizing the 
magnetic field changes, for a current density vector constrained to belong to 
some bounded set (J A). It is shown that the selection of A determines 
the specific critical state model in use. Meaningful choices of A are discussed 
in view of the related physical mechanisms. 

PACS numbers: 4I.20.Gz,74.60.Jg, 74.60. Ge, 02.30.Xx 



1. INTRODUCTION 

Macroscopic magnetic experiments, which are sensitive to thermody- 
namic averages over the sample's volume, have been a fundamental tool in 
the research of type-II superconductors. In particular, one can obtain infor- 
mation on the interaction between flux tubes (vortices) and the underlyinff 
pinning structure. Within the framework of Bean's critical state modellil 
(CSM) the experimental data can be processed with remarkable ease. Thus, 
elementary mathematical manipulations are usually enough for getting im- 
portant material parameters as the critical current density Jc. 

Being a phenomenological approach, the CSM is supported by meso- 
scopic theoriesp which deal with the complex vortex interactions. The phys- 
ical basis is that a maximum gradient in the density of vortices is allowed by 
the underlying pinning structure. This gradient is nothing but the critical 
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value of the coarse-grained current density. When more and more vortices 
are pushed inwards by the external field pressure, avalanche processes are in- 
duced as Jc is exceeded, which relax to new metastable structures with max- 
imum current density. Overdamped relaxation is typically assumed, which 
relates to the fact that external drive variations occur slowly as compared 
to the superconductor's diffusion time. 

From the macroscopic point of view, this course of events may be mod- 
elled by the use of the Maxwell equations and a material law in the form of 
a vertical E{J) characteristic (see Fig|l|). Any dc current density may flow 
within the sample in the absence of electric field for the range — Jc < J < >/c- 
Corresponding to a very high level of dissipation for \J\ > Jc, magnetic field 
changes are introduced by electric field values along the vertical line. Even- 
tually, a new dc state is reached when E comes to zero again, and we have 
J = lb Jc in those regions of the sample affected by changes. 

For one dimensional problems with parallel vortices, the actual behavior 
of E may be skipped over, and many works on the critical state do not even 
mention this quantity. In fact, the solution of the CSM consists of choosing 
among one of the possible states for the current density {— Jc, 0, Jc} accord- 
ing to the rule that any electromotive force due to external field variations 
induces the maximum current density flow. Nevertheless, as soon as mul- 
tidimensional situations are considered (nonparallel vortices), the problem 
becomes much more difficult because an infinite number of possible states 
for the critical current density Jc are allowed. Several approximations have 
been introduced over the last decades for dealing with such a situation. 

In the early 70's the consideration of simultaneous application of longi- 
tudinal magnetic field and transport current along type-II wires led to the 
concept of flux compression between mutually perpendicular vortices. This 
approach allowed to explain the so-called paramagnetic effect related to the 
increase of magnetic flux towards the center of the sample and basically 
neglects the flux cutting events between nonparallel vortices.3 

However, in the following decade, flux cutting was recognized to be 
the essential mechanism for a reasonable understanding of rotation field 
experiment^li'i and gave place to the concept of double critical state modefl 
(DCSM). In this approach, the so-called parallel critical current density J^\\ 
and perpendicular critical current density Jc^ were introduced. J^n relates 
to the flux cutting threshold, whereas Jc_l stands for the conventional de- 
pinning current density. 

Recently, a more sophisticated approach was introduced, the so-called 
two-velocity hydrodynamic modelu Essentially, this theory incorporates the 
flux pinning and cutting phenomena within the framework of flux line lattices 
(FLLs) which consist of two vortex subsystems, moving at different speeds. 
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This physical description relates to the mathematical need of a discontinuity 
in the velocity field in order to allow for crossing events. 

Finite element based models are also available,El which allow to compute 
magnetization properties of many non-ideal geometries of crucial importance 
in the field of applied superconductivity. These are typically in the form of 
variational statements and, to our knowledge have been exploited for scalar 
E{J) relations (i.e.: ii^ || J is assumed). 

A variational approximation to the critical state was also suggested in 
Ref.p!0|, by a guessed minimum flux change criterion. The authors of that 
work argued that, in certain conditions {E || J), successive field penetration 
profiles in the magnetization process could be obtained by minimizing the 
squared flux density variation for constant current density modulus. 

In the theory presented in this paper we propose a generalization, which 
incorporates a vector vertical E{J) law, but in the manner of Bean's ap- 
proach still avoids the explicit use of induced electric fields. A principle is 
provided for determining the critical current density vector Jc within the 
sample, subsequent to the breaking of equilibrium by the external drive. We 
keep the main physical ingredient that infinite dissipation occurs if some 
critical value is exceeded. As this value is a vector quantity a new mathe- 
matical framework is required, so as to find the actual critical state in which 
the system settles subsequent to a specific manipulation. 



2. MACROSCOPIC FIELD EQUATIONS 

In this section we will define the macroscopic fields involved in the CSM 
and their relation to measurable quantities. Our work has been developed 
within the assumption that the sample's response is dominated by the ir- 
reversible contribution. This means that the equilibrium response of the 
FLL in the absence of pinning centers is neglected. For global magnetiza- 
tion measurements in which the integrated magnetic moment of the sample 
is measured, this is a very good approximation, unless for weak pinning 



specimens (see Ref.ll and references therein). Then, the coarse-grained 



electrodynamics if formulated in terms of 

(i) The flux density B{x) within the sample is the average of the micro- 
scopic field intensity h. 

(ii) The magnetic field H{x), in the absence of equilibrium magnetization, 
is linearly connected to B hy B = ijlqH. 



(iii) The averaged current density may be calculated from J = W x H . 
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Fig. 1. Restriction sets for the current density (J G A) 

(iv) The electric field can be evaluated by V x ^ = —nodtH. 

Finally, the connection to observable quantities is done by characterizing the 
external field sources and the sample's response. 

(v) On neglecting finite size effects, the magnetic source {Hs) enters as a 
boundary condition for the flux density at the surface of the sample. 
The tangential component is continuous BT{surface) = fioHs^r- 

(vi) Then, the measured magnetic moment of the sample per unit volume 
is M = {B)/fio -Hs = {H) - Hs. 

For further use we notice that, in the case of one independent variable. 
Ampere's law allows the following representation 

-> -> dH da ^ 

j = ^xH = —a-H—H = J±a + JiiH (1) 
dx dx " 

where a is the angle between the magnetic field and the Z axis, and H and 
a are unit vectors respectively directed along H or perpendicular. 



3. OPTIMAL CONTROL THEORY FOR THE CRITICAL 

STATE 

In this section we show that Bean's ansatz for parallel FLL's allows 
a geometrical interpretation, which one can straightforwardly generalize for 
complex systems with nonparallel vortices. Recall that penetrating flux pro- 
files in the slab geometry may be obtained by solving the equation 
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[n5t,H(x)] 

Fig. 2. Quasistationary evolution approximation 

where the sign selection is made according to the rule: any electromotive 
force due to external field variations induces the maximum current density 
flow. 

In simple terms, one is just using Lenz's law. Thus, the macroscopic 
electrodynamics of a hard type-II superconductor may be viewed as follows. 
Owing to physical limitations imposed by the material properties, the cur- 
rent density could take any value within the interval A = [— Jc, Jc], but as 
electromagnetic induction operates all along the process, one gets J = ±Jc, 
according to the most effective way of shielding field variations. From our 
point of view, the physical limitations within the superconductor define a 
region in J-space where the current density could settle (J G A). Lenz's law 
establishes that only the boundaries will be occupied (or the point J = for 
regions where electric fields have not reached). 

For higher dimensional systems, i.e.: vortices are no longer straight par- 
allel lines or J is truly a vector [see Eq.(|l|)], the allowed region for the current 
density will be a subset of J-space. In this case, Lenz's law is insufficient 
for determining the points where the current density lies along the magne- 
tization process. Faraday's law in full form must be appealed. In brief, our 
critical state theory includes (i) a material law of the kind J € A and (ii) a 
variational statement of Faraday's law. Below, it will be shown that calling 
on (ii) , the solution of the problem J * only takes values on the boundary of 
A (or 0). This has been sketched in Fig.|l[ Thus, the essence of the critical 
state in superconductors, is defined by the condition 

J* G OA (or 0) . (2) 

The determination of the specific values J*{x) for any process is described 
in the following subsection. A brief discussion on the physical basis of the 
mathematical statement is also given. 
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3.1. General statement 



As it was shown in previous workH^fiJ, the evolutionary magnetic pro- 
cesses in hard type-II superconductors may be described by quite general 
principles of linear irreversible thermodynamics. Let us consider a supercon- 
ductor in some equilibrium state within the experimental range of interest. 
This means that the system of vortices is at rest by compensation of mu- 
tual interactions, pinning forces and the external field pressure. From the 
macroscopic point of view, the equilibrium states will be characterized by a 
magnetic field or current density pattern within the sample. For instance, 
one can define it by a vector function Jc . (x) . Consider a small perturbation 
that drives the system out of equilibrium. Then, a non-equilibrium flow 
and linearly related internal forces will arise. Specifically, the critical cur- 
rent density has been exceeded by some amount Jj — Jc-, linearly related 
to associated electric field components Ej. Stationary states in which the 
external drive keeps time independent values of E can be straightforwardly 
solved by the principle of minimum global entropy production, introduced 
by Prigogine in the context of linear irreversible thermodynamicsB. In this 
case, one can write it as 

(e-{J-J,) ^ min, (3) 
Jn 

where i7 denotes the volume of the sample. The reader can check that in 
the particular case Jc = 0, this principle leads to the well known current 
distribution for the stationary regime in normal metals. El 

However, we are interested in the time evolution between equilibrium 
states induced by successive magnetic field changes. One can achieve an 
appropriate description just by adding up a term to Eq.(P) which accounts 
for the magnetic field inertia, in the so-called quasistationary approxima- 
tion. The physical idea is schematically shown in Fig.|2[ We assume a time- 
discretized picture in which n5t denotes the n-th time layer after n equal time 
steps with increment 5t. Then, for small enough values of 6t, one can also 
suppose that changes in the magnetic field configuration H^{x) are small and 
substitute the continuum trajectory of the system by a polygonal approxi- 
mation in which (i) H^[x) changes at fixed time to H^j^i{x) [inertia term) 
and (ii) time is incremented for Ei = Pij{Jj — Jcj) = Pij[C^ x — Jcj] 

(stationary entropy production term). Then, the configuration H^-^i{x) may 
be obtained from the previous one by the minimization rule 



fio \H^+i - + St E ■ (J - J,) ^ min. (4) 
Jn Jn 

The rigorous proof that the previous principle is equivalent to Faraday's law 
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in discretized form plus a constitutive relation of the kind Ei = Piji Jj — Jcj) 
is straightforward by following the calculations in Ref.|l^. 

The final form of our principle, which connects to the geometrical in- 
terpretation described in the previous paragraph is attained by considering 
vertical E{J) laws, i.e.: infinite flux flow resistivity. Then, in order to avoid 
adding a divergent positive quantity, the entropy production term is replaced 
by the constraint J G A. We end up with the principle 

In a type- II superconducting sample with an initial field profile H^{x), 
and under a small change of the external drive, the new profile IIn+i{x) 
minimizes the functional 

C[H,,+i{x)] = I [ \Hn+i - HnW (5) 

with the boundary conditions imposed by the external source, and the 
constraint V x H^+i G A. 

In general, A is a bounded region where the dc current density vector 
should remain. As it will be shown in the next subsection, the selection of 
this region corresponds to using a given critical state model. 

The kind of mathematic al sta tement posed in Eq.(^) may be solved by 
several methods. The authorslSOJ have preferred to use the so-called Opti- 
mal Control (OC) theory,E3 a powerful generalization of variational calculus, 
which allows to deal with a wide class of constrained problems. For the 
readers' benefit we have included a brief survey of this mathematical tool 
as an appendix of this article. A reminder of the application to the critical 
state problem is given below. 

For the slab geometry, with applied magnetic field parallel to the faces, 
the current density restriction may be written 

^^ = /(i?„+i,n,x), /gAx. (6) 

This is called control system. Hereafter, we take the X axis perpendicular 
to the slab faces and the origin of coordinates at the midplane. By con- 
struction, the vector / = {0,fy,fz) is orthogonal to the physical variable 
J, fy = Jz, fz = —Jy Notice that, owing to the applied rotation, we use 
the notation A_|_ for the allowed control set. Note also that the function f 
displays dependence oiiJJie local magnetic field (this allows to include the 
usual models Jc(|-^|),t3't3 and possible anisotropy), on the position (for po- 
tentially inhomogeneous materials) and on some independent coordinates H, 
the control variables in OC terminology, parametrizing the region A^. 



A. Badia and C. Lopez 





A 


P 














f 







Fig. 3. Maximum projection rule (=^ J* € dA) 



Next, we define a Hamiltonian density, containing the Lagrangian den- 
sity to be minimized 

n{Hn+l,U,P,x) = p- f- ^\Hn+l - Hnf' ■ (7) 

Denoting by H*_^,l{x), p*{x) and u*{x) the optimal solution functions (i.e. 
minimizing C and satisfying the control system), the OC equations are 

^^ = ^ = /(^„+i,n ,x), (8) 
the adjoint equations for the momenta 

ax dHn+i dHn+i 
and the algebraic condition of maximality 

n{H,u*,p,x) >n{H,u,p,x) y f{H,u,x) £ A± . (10) 

For the class of Hamiltonian Tl described above, the algebraic condition 
of maximality is fulfilled for a vector / with maximum projection over the 
momentum p. As it apparent in Fig.^ and in the various examples of the 
following subsection, this rule is nothing but the result announced in Eq.(|2|): 
J* G dA. The actual solution J *{x) can be determined by the Hamiltonian 
equations displayed above. 

As regards the boundary conditions that must be used for obtaining 
the solution profiles H*^^{x) from Eqs.(|8|) and (^), several considerations 
must be made. Firstly, in the absence of demagnetizing effects, H*^^ is 
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Fig. 4. Critical State Models arising from the selection of the region A 



determined on the faces of the slab by continuity of the external applied 
field. For the remaining boundary conditions, two typical situations appear. 
In the first case, the modified penetrating profile equals the former 

profile Hn{x) before reaching the centre of the slab. This gives place to an 
(unknown in advance) point x* such that H*_^_i{x) = Hn{x) V x < x* . 
The algebraic free boundary condition TC{x*) = applies in this case (see 
[Appendix A. ), and allows to determine x*. In the second case, the new 
profile never meets the former one, and the value of H*_^i{0) is unknown; 
then, the corresponding transversality condition p{0) = (see also Append!^ 
A.) completes then the number of required boundary conditions. 



3.2. Particular models 

Next, we analyze the models which arise for several selections of the 
restriction region in Eq.(^. FigJ^ depicts such regions and the relation of 
the functions / with the associated momenta in the theory. In this section 
we will focus on the formal aspects of the solutions. Later on, the relation 
to the experimental facts will be established. 
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3.2.1. Isotropic model 

The simplest restriction set which one can imagine for two dimensional 
systems corresponds to choosing A as a circle. Following the notation intro- 
duced in Eq.(P), the control system takes the form 



dx 



f{Hn+i,u,x) = u . (11) 



In this case, n is a vector within the unit disk D. For simplicity, dimensionless 
units have been used: x is given in units of the slab half-thickness a, J in 
units of Jc, and the magnetic field H in units of JcC. Also for simplicity, 
here and in what follows we take Jc as a constant. 

The maximality condition (maximum projection of / over p ) leads to 
the solution /* = u* = p* /p*. Note that this means that /* is a unit 
vector always parallel to the momentum as indicated in Fig.^. Eventually, 
the penetration profiles H^_^^ j(a;) can be obtained from the system 



dH, 



^ = (12) 
dx p* 

dp* 



dx 



^n+l,i - -f^n.i (13) 



on supplying the appropriate boundary condictions. 

The most relevant physical property related to this model is that the 
maximum allowed current density modulus Jc is carried within those regions 
which have been affected by the perturbation ( J * G dA in our geometrical 
language). Then, when a considerable amount of current is required along 
a given direction (f.i.: J^) for shielding a specific field variation (Hy in this 
case), the other component of J is depressed {Jy in this case). This reduces 
the shielding capability in the orthogonal field direction (Hz in this case). 
As a consequence, one can effectively increase the penetration of a given 
magnetic field component, just by manipulating the nerpendicular one. This 
phenomenon has been named magnetization collapse^ and has been widely 
illustrated in the framework of our theory in Refjl^. 



3.2.2. Double critical state model 

Below, we show that the DCSM equations can be readily obtained 
within the framework of our theory. A specific selection of the restriction 
region A, together with the variational principle, render the DCSM as a 
particular case of our formulation. 
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Recall that the flux depinning and cutting thresholds are defined by the 
critical state conditionsB 

,dH^+i dan+i. , . 

\^^\<Jcl. , l^l<^c||. (14) 

Here, a{x) is used for the angle between the local magnetic field vector and 
a fixed axis. 

On making use of the possiblity of expressing the Lagrangian in conve- 
nient generalized coordinates, the OC functional may be written 

C[Hn+l{x)] = \ I -ffn+1 - 2HnHn+l COS (On+l - On) • (15) 

2 Jn 



The modulus-angle representation enables to write the statement of Eq.(14) 
as 



J± = f± ; —7— = k\\ = /|| 



dx dx 
with (/||,/_l) a vector belonging to a rectangular region of size 2kc\\ x 2Jc± 
(see Fig.^). Notice that, in this case, the region A itself is depicted. We also 
recall that, on using the appropriate dimensionless units A may become a 
square. Thus, if one defines A = l/A^^n as the length scale and Jc_lA as the 
magnetic field unit, / becomes 

/± = Uh ; \uh\ < 1 

f\\ = Ua ; \Ua\ < 1 . 

For simplicity, this convention will be used below. Then, the associated 
Hamiltonian is 

W = PhUh + PaUa - ^ [Hn+l " 2H^H^+i COS (ttn+l " "n)] , (16) 

and the maximality condition is again fulfilled by J * G 9A. However, as 
one can also notice in FigJ^, the rectangular shape determines new features 
in the critical current density vector behavior. If the momentum vector p 
has non-vanishing components one gets /* = u* = [sgn(pQ), sgn(p/j)], i.e.: 
/ * lies in a corner of the rectangle. On the other hand, if either pa or p^ 
vanish, / • p maximality occurs for any vector leaning somewhere on a given 
side of the rectangle. Then, / * must be determined calling on some addi- 
tional condition. As we have indicated in Fig.^ the previous mathematical 
properties relate to the physical concept of CT, T and C zones introduced 
by Clem and Perez-Gonzalea3 . Thus, starting with Eq.([l6|) one gets 



(17) 
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sgn(po) (18) 

i/n+l - -H'n COS (On+i - On) (19) 

i7n+ii^nsin(an+i - On)- (20) 

The solutions of these equations may be classified in terms of the momentum 
variable behavior 

(i) CTzone (p„/0,p;,/0) 

The field penetration is given by 

Fn+i = sgn{ph)x + Hj,+i{0) 
On+i = sgn(p„)x + an+i(0) . 

(n) T zone {pa = 0,Ph7^ 0) 
We have 

-f^n+i = sgn{ph)x + Hn+i(0) 
and an+i{x) is to be determined by the condition dpa/dx = 0. Then 

H,,+i{x)Hnix) sm[an+iix) - a^ix)] = . 

(in) C zone (pa 7^0,ph = 0) 
We have 

On+l = Sgn{pa)x + an+l(0) 

and Hn+i{x) is to be determined by the condition dph/dx = 0. Then 
Hj,+i{x) = Hn{x) cos [an+i{x) - Onix)] . 

We notice that the modulus and angle profiles display maximum (criti- 
cal) slope for the so-called CT zone, corresponding to the fact that both the 
depinning and cutting thresholds have been exceeded in that region. Thus, 
fiux Cutting as well as flux Transport have occured. Within the so-called 
T and C zones one gets maximum slope for one of the two components of 
the field and subcritical behavior for the other. In physical terms, the highly 
resistive state has only been triggered for a particular direction in J-space. 

Finally, we want to mention that the full equivalence of our interpre- 
tation to the standard DCSM formulation has been tested.113 In fact, it has 
been shown that the latter is nothing but the continuum limit of our time- 
discretized approach. 



dx 

dph 
dx 

dpa 

dx 



The Critical State in Superconductors with Cross-Flow 




Fig. 5. Rotation angle and current density modulus within the elliptic model, 
a is given in radians and J in dimensionless units 



3.2.3. Elliptic model 

It is well known that type-II superconductors may display anisotropic 
physical properties, either related to their crystal structure or artificially 
induced by a number of processes as thermal or mechanical treatment, irra- 
diation, etc. Thus, we are led to consider the influence of this phenomenon 
on the critical state. The simplest assumption which one can make is that the 
modification of the pinning strength along perpendicular directions within 
the sample can be described by an elliptic region A (see Fig.^). At least, for 
moderate anisotropy, this may be considered as a first order perturbation 
of the isotropic model, which provides a good phenomenological theory for 
many experimental facts. 

From the mathematical point of view, the OC equations are also a 
straightforward modification of the isotropic expressions. Thus, the control 
system becomes 



C^^n+l,y 

dx 
dx 



fy = lUy (21) 
f.=u,, (22) 



with 7 the anisotropy parameter and H a vector within the unit disk. This 
leads to the associated Hamiltonian 

n = Py'-fUy + PzUz - ^(^n+l - -^n)^ • 



In this case, the maximality condition produces new features in the optimum 
solution. As one can notice in Fig.^, maximum projection does no longer 
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mean / || p as in the isotropic case. In general, these vectors will be at 
an angle < (3 < vr/2, which can be close to 7r/2 for 7 far from unity. 
Specifically, we get 

J I ? 

and, thus, the Hamiltonian equations read 



(23) 
(24) 



dx 
d^ 

dx 



Hn+l,y - Hn,y (25) 
Hn+l,z - ^n,z • (26) 



Notice that this system includes the isotropic model as a particular case 
(7 = 1) and that excitation fields along the principal axes (Y and Z in this 
case) produce one- dimensional critical state models with unequal critical 
currents as expected. 

For illustration of the field penetration properties in this anisotropic 
model, we include Figjs], in which some aspects of a field rotation experiment 
are depicted. We have simulated a field consumption process for a field 
cooled type-II slab with nonmagnetic initial state {Hq = constant z) and 
7 = 0.5. Notice that, as the magnetic field vector rotates within the sample 
the current density modulus J{x) varies from 0.5 to 1. 



3.2.4- Pseudo-isotropic model 

In this part, we develop another approximation which is also suited 
for dealing with anisotropic systems. Within the so-called pseudo-isotropic 
model, we assume that the restriction region A is a circle, but with vari- 
able radius according to the magnetic field orientation with respect to some 
particular axis within the sample (see Fig.^. This can be written in the 
form 

|J| < Jc/(a), 

with Jc a constant and /(a) some function of the angle between the local 
field and the Z axis. In order to allow comparison with the elliptic model. 




we introduce the control system 



C?-ffn+l,y _ 7-^n+l,y + ^n+l,z 

dx ~ Hl^^ 



= fy = fi<^)^y (27) 



— ^ — - = ^72 u^= fz= f{a)uz ■ (28) 

ax -"n+l 

Once more, u is a vector within the unit disk. It is apparent that for the 
particular case of parallel FLL's with magnetic field along principal axes, one 
gets the same one- dimensional approximations which arise from Eqs.(21 ) and 
(22). For instance, = ^ fy = 'yuy, fz = 0. 

The minimization of the cost functional C is now attained by maximizing 
the Hamiltonian 

ri - 772 P ■ u - -{Hn+l - J^nj ■ 

This leads to the condition /* = f{(x)p*/p* and one gets the Hamiltonian 
equations 



^-^n+l,y _ 7-^n+l,y + 

r*2 



(29) 



dx p* H*^^ 

(30) 



dH*+l,. ^ pl 7i^nA,y + gnA 

dx p* H* ^ 



'n+l 

dP*y _ rr* rr (7 - l)g^+l.ygn*|l,z . 

— - H^^,^^-H^,,-2p (31) 
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dpi _ rr* rr n^* ~ 7)-^n+l,z-^n+l,y , . 

- ^n4-l,z - ^n,. - 2p . (32) 

As in the previous model, we illustrate some properties of the magnetic 
field penetration. Fig.^ depicts the current density variation as the local 
field rotates within the sample for a nonmagnetic initial state experiment as 
in Fig.p. Notice that the current density modulus also bounces between the 
values 0.5 and 1, but with a clearly different structure. Although both models 
share the same one- dimensional limits, the actual current distribution J{x) is 
influenced by different factors. Thus, within the elliptic model, the shielding 
current value is determined by the external fleld variation. For instance, for 
the initial rotation steps one gets J ~ 0.5 because the induced current is 
mainly formed by Jy. On the contrary, the magnetic response within the 
pseudo-isotropic approach is determined by the angle a itself. The actual 
orientation of J is not important. Thus, the initial rotation steps are shielded 
by J ~ 1 because H is still basically directed along Z axis. 



4. COMPARISON TO EXPERIMENT 

In this section we will discuss the application of our theory for the 
analysis of some relevant experiments related to the investigation of cross- 
flow effects. According to the boundary conditions induced by the external 
drive, they may be classified as: (i) rotation experiments and (ii) crossed- 
field experiments. In the first case, the applied magnetic field Hg rotates 
at some constant and low enough frequency ojq, so as to neglect possible 
relaxation effects. The field modulus Hs remains constant. In the second 
case, the orthogonal components of Hs are cycled for some definite process. 
For instance, one of them is raised and lowered while the other one is kept 
constant. 



4.1. Rotation experiments 

The main features of rotating field experimentsSHi have been well de- 
scribed within the original formulation of the DCSM.@ For instance, in the 
so-called nonmagnetic initial state, the magnitude of the fiux density within 
the sample is decreased as if vortices were somehow expelled from the sam- 
ple. Indeed, Clem and Gonzalez clarified that such effective expulsion may 
be naturally understood as a fiux cutting phenomenon between adjacent ro- 
tating vortices. Within the DCSM, the following picture was given. The flux 
density develops a diamagnetic profile toward the sample's midplane, until 




it becomes zero for some distance xq (or until the midplane is reached). In 
the case that xq > exists, the field modulus profile reaches a stationary 
V-shape as soon as xq is defined. 

We showed that the behavior described above mav also be obtained 
in the framework of our theory for the isotropic model£3 It was remarked 
that the decoupling point xq is sharply defined by the process itself, and 
that further rotation leaves the field modulus unchanged. This behavior 
has been emphasized here (see Fig.^). We have plotted a number of profiles 
subsequent to the issue of xq. Notice that Hy{x) and Hz{x) are frozen for x < 
Xq. Notice also that, though these components evolve for x > xq, one gets 
a frozen modulus profile and nearly linear penetration of the rotation angle 
a{x). Physically, the process may be described as the separation of vortices 
in two groups, one which rigidly pins to the sample (inner vortices), and 
another one which frictionally rotates relative to the sample while keeping a 
constant density profile (outer vortices). 
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Fig. 8. Simulated field penetration profiles for an anisotropic sample under 
external rotating magnetic field 

Finally, we will show that, as it was already guessed in Ref.^, the ap- 
pearance of higher frequency oscillations (2^0) in the magnetic moment com- 
ponents can be ascribed to anisotropy in the current flow. Fig.|8| depicts our 
simulations for the field penetration profiles in a rotation experiment as 
calculated by the pseudo-isotropic model defined in the previous section. 
We want to notice the appearance of a V-shape modulus profile. However, 
by contrast to the isotropic case, this structure keeps no longer station- 
ary. Further rotation produces a complex two minima structure. A second 
turnabout (27r < uq < An) defines a new V structure and the process goes 
on in a periodic fashion. As regards the angle penetration profile a{x), a 
multiple step structure is observed. The definition of decoupling points is 
accompanied by jumps in the rotation angle. Subsequent steps produce a 
negative angle for inner points, which is interpreted as a counter rotation of 
flux lines. The double frequency signal, which was reported in the early 80 's 
by Boyer and co-workers is straightforwardly related to the physical mecha- 
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Fig. 9. Simulated magnetization components for an anisotropic sample under 
external rotating magnetic field 



nisms described above. In Fig.^ we have plotted the averaged magnetization 
components both longitudinal (Ml) and transverse (Mt) to the applied field 
Hs- The plot collects the evolution Ml^t{os) with respect to the applied 
rotation angle for three different values of the field modulus (0.1 ,0.25 ,1 
in our dimensionless units). Note that, for the lowest field value, M^t^ 
basically display a stationary harmonic behavior at the sample's rotation 
frequency ujq. As increases, one can readily observe a frequency mix- 
ing phenomenon. Eventually, for high enough fields, M becomes a rotating 
vector at frequency Iluq. 

The rotational scenario of field cooled anisotropic samples with non- 
magnetic initial state is established as follows. Two groups of vortices exist: 
a rigid core well within the sample and a region of fiux tubes below the 
surface, which undergo complex transport and cutting phenomena as ro- 
tation proceeds. The rigid core contributes as a harmonically oscillating 
magnetic moment when analyzed in a reference frame at rest with respect 
to the applied field. On the other hand, the active region of vortices close 
to the surface contributes with an essentially double frequency signal. For 
low enough fields {Hs = 0.25) the rigid core nearly spans over the whole 
sample and the response is linear. For high enough fields {Hs = 1) the rigid 
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Fig. 10. Simulated magnetization processes for crossed field experiments in 
the configurations used by Park and Fisher 

core has disappeared and the response is at 2ujq. Intermediate values hold a 
frequency mixing, corresponding to the superposition of both effects. 



4.2. Crossed field experiments 

For illustration of the characteristic phenomena in the crossed field con- 



figuration, we will discuss the results within Refs.1^,1^. Hereafter, they will 
be named Park's and Fisher's experiment respectively. In both cases, the 
superconductor is firstly subjected to a magnetic field along a given direction 
(f.i.: Hzs)- Subsequently, this component is fixed and the orthogonal one 
(Hys) is cycled while recording the magnetic moment. 

Park and co-workers measured My and for high-Tc crystals. Their 
results show that vanishes as Hys is cycled. Additionally, a more or 
less conventional loop My{Hys) is observed. However, the loop fails to close 
after a field cycle has been completed. Both features are reproduced by the 
application of our theory within the isotropic model (see Figjl0|) . 

Finally, the symmetric supression of the magnetic moment by cy- 
cling Hys reported by Fisher may also be explained within the isotropic 
hypothesis. In Fig.|l^ we display our simulations of the diminishing behav- 
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ior of Mz, both for the para- and diamagnetic states. In the first case, Hzs 
was cycled up to a positive magnetic moment state, and then fixed. In the 
second case, HzS was raised up to a definite negative magnetic moment state. 

The underlying physical mechanisms responsible for the described be- 
havior have been extensively discussed in Ref.l3. Basically, owing to the 
rectriction \J\ < Jc, the shielding capability related to is decreased when 
a considerable amount of current is required for minimizing /(//n+i,y— -f^n,y)^) 
which is the dominant process in the configurations discussed. 



5. CONCLUDING REMARKS 

The fundamental subject of our work is that a number of observables 
related to experiments with crossing vortices in hard superconductors may 
be well described in terms of critical state principles. In particular, we 
have shown that simulated rotation and crossed field experiments both for 
isotropic and anisotropic systems nicely reproduce the actual measurements. 

It is shown that Bean's model for parallel vortices can be thought as a 
particular case of a general critical state theory. The main ingredients of the 
theory, which is given in the form of a variational principle, are the magnetic 
field inertia (Faraday's law) and a very high level of entropy production for 
moving vortices when the threshold for static configurations is exceeded. 
The arising constrained minimization problem has been solved by optimal 
controlS'ij techniques, which nicely fit the physical statement. However, we 
want to notice that other methods can be applied to solve the problem posed 
in Eq.(^). In particular, we suggets the use of the so-called direct method^ 
in which a family of appropriate base functions for the solution generates a 
problem in matrix form as in standard finite-element theories. 

Although many relevant experimental facts can be explained within our 
theory, several approximations have been used in order to keep the mathe- 
matical intricacy to the lowest degree. An improved quantitative description 
would require the consideration of finite size effects, the influence of equilib- 
rium magnetization, and the finite value of the flux-flow resistivity. These 
are merely technical aspects which can be incorporated by using OC for par- 
tial differential equations, a definite B{H) relation and a dissipation term in 
the variational quantity. 
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Appendix A. A BRIEF SURVEY ON OPTIMAL CONTROL 

Control theory is a branch of dynamical systems theory in which the 
usual ingredients are a family of evolution differential equations of the system 
containing the time, some phase space coordinates representing the states of 
the system, and extra control variables modelling some external action ap- 
plied over it. The applications of the theory have been mainly developed in 
the field of Engineering. Within this realm, the controls are the representa- 
tion of some manual or automatized action and, thus, affected by limitations. 
In mathematical terms, the controls belong to a bounded set. 

Self-organized critical systems, which are a flourishing research area in 
Physics, seem to be a natural candidate for the application of this mathe- 
matical tool. In this case, the control variables and limitations are related to 
physical interactions. In particular, the magnetic properties of hard type-II 
superconductors rely on the interaction between flux tubes and a maximum 
force from the underlying pinning structure. 

In order to develop the optimal control theory in a brief but self- 
contained presentation, one should begin by recalling some concepts on the 
classical variational calculus (see Ref.^ for a comprehensive and detailed ex- 
position of the topic). Civen a Lagrangian function of position and velocity 
L{x,dx/dt), the problem of minimizing the action integral jj^ L{x,dx/dt)dt 
for a curve x{t) with fixed endpoints x{a) = Xa, x{b) = X{, generates, by 
applying a first order variation 6x{t) of the curve and after integration by 
parts, the well-known Euler-Lagrange equations 



Above, the boundary conditions which imply the conditions 6x{a) = 6x{b) = 
have been invoked in the integration by parts. For other boundary condi- 
tions, these terms do not vanish authomatically, and some extra transversal- 
ity conditions appear. For example, if the final point is not fixed but belongs 
to some hypersurface S C IR", then its first order variation 5x{b) is tangent 
to the hypersurface, and therefore the extra condition is 



The Hamiltonian version of the classical variational calculus, so impor- 
tant in the process of quantization of the fundamental physical theories, and 
also in the study of many dynamical systems, will be obtained below using 
a kind of trick that will be useful later on. Let us consider that coordi- 
nates V are independent of coordinates x, and afterwards let us impose the 
constraint v{t) = dx/dt. In this case, on using the method of Lagrange 
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multipliers for constrained minimization problems, we must use the action 
integral jj^ poL{x,v) + p ■ {dx/dt — v)dt, where p are the Lagrange multiplier 
variables, and po is a positive constant for minimizing problems. The Euler- 
Lagrange system for this extended Lagrangian comprises three families of 
equations, one for each base variable {x,v,p). We get 

dpi dL 
dt dx^ 



for first order variations in x, 



= ^ 
dt 



for first order variations in p, and 

= — -Pi 

for first order variation in v. 

The last equations define the Legendre transformation between velocity 
{x,v) and phase {x,p) spaces. For most cases in the context of classical 
mechanics, this map is an isomorphism and can be inverted to solve for v 
as a function of x and p. This relates straightforwardly to the condition 
deiip^L/dv'^dv^) / 0. Then, it is a simple exercise to check that the other 
two families of equations can be written in the form 

dpi _ _dH dx' _dH 
dt dx^ ' dt dpi 

where H{x,p) = p ■ v{x,p) — L(x, v{x,p)) is the usual Hamiltonian. 

For some pathological cases, the Lagrangian function is not regular and 
the Legendre map is not invertible. However, a Hamiltonian theory can still 
be constructed, but the true Hamiltonian equations live in a mixed velocity- 
phase space, with independent {x,v,p) coordinates. In such problems, one 
of the families of equations (those determining v ) is not differential but al- 
gebraic. One has a so-called DAE system (Differential Algebraic Equations) 

dpi _ _dH dx' _dH - 

dt dx^ ' dt dpi ' dv'^ 

with H{x, v,p) = p- V — L{x,v). 

Let us now introduce a control system for the velocities 

dx - ^ 
^ = /(x,«). 
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Here, the coordinates u are the so-called controls, in general, taking values 
within some subset A C IR'^ . These extra coordinates may be modelling some 
external action on the system or also some intrinsic limitation as it is the case 
of the critical current density in superconductors. Let us also assume that we 
want to minimize a functional J L{x,u)dt. We notice that the control system 
may be thought as a parametric description of some subset of the velocity 
space, and therefore as a constraint for the minimization problem. This 
constraint directly relates to the singular behavior of the matrix d'^L/dv'^dv^ 
and the method for pathological systems must be applied. One must consider 
the Hamiltonian equations for H{x,u,p) = p ■ f{x,u) — L{x,u), i.e. 

dx* dH : ^ 

dpi _ dH _dL ^ df 
dt dx^ dx^ ^ dx^ 
Q _ dH _ ^ df dL 

Note that the last equation is a first order necessary condition of minimality 
relating only to stationary internal points of the region A. For general 
problems with bounded control spaces, a stronger condition for determining 
the solution u* isrequired. This is given by the so-called maximum principle 
of Pontryagiir^^n, which we state without proof here 



u* B H(x, u* ,p) = max H(x, u,p) ■ 



As far as this condition can be explicitly solved, we get u*{x,p), and after 
back substitution into the other equations we get a system of first order 
ordinary differential equations. As regards the boundary conditions, they 
can be simple as in the case of fixed endpoints for x, or they can include 
transversality conditions, as pointed before. For example, when the final set 
for X is an hypersurface of IR", the transversality condition means that f) 
is perpendicular to the hypersurface (in particular, free final point x means 
that the hypersurface is IR" itself and implies a vanishing final momentum). 
Sometimes the problem does not have a defined final point 6; then, the 
associated transversality condition is the vanishing of the Hamiltonian at 
this unknown point .c3 

Just for simplicity, the previous exposition has been done within the 
scenario of particle dynamics. Of course, it can be extended to continuous 
systems by using the fields as dynamical variables and defining the corre- 
sponding Lagrangian and Hamiltonian densities. 
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